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Abstract. Self-gravitating systems such as elliptical galaxies appear to have a constant integrated specific entropy 
and obey a scaling law relating their potential energy to their mass. These properties can be interpreted as due 
to the physical processes involved in the formation and evolution of these structures. Dark matter halos obtained 
through numerical simulations have also been found to obey a scaling law relating their potential energy to their 
mass with the same slope as for ellipticals, and very close to the expected value predicted by theory. Since the X- 
ray gas in clusters is weakly dissipative, we test here the hypothesis that it verifies similar properties. Comparable 
properties for the dark matter component are also investigated. 

With this aim, we have analyzed ROSAT-PSPC images of 24 clusters, and fit a Sersic law to their X-ray surface 
brightness profiles. We found that: 1) the Sersic law parameters (intensity, shape and scale) describing the X-ray 
gas emission are correlated two by two, with a strong correlation between the shape and scale parameters; 2) the 
Oh' hot gas in all these clusters roughly has the same integrated specific entropy, although a second order correlation 

between this integrated specific entropy and both the gas mass and the dynamical mass is observed; 3) a scaling 
j.^ ■ law links the cluster potential energy to its total mass, with the same slope as that derived for elliptical galaxies 

and for dark matter halo simulations. Comparable relations are obtained for the dark matter component. All 
these correlations are probably the consequence of the formation and evolution processes undergone by clusters 
of galaxies. 
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1. Introduction cussed by Gunn & Gott (1972), and self-similar solutions 

for dark matter halos and gas were studied in numeri- 

Clusters of galaxies are known to be the largest gravita- ^ simulations carried out e g by Fmmoie & Qoldreich 

tionally bound objects in the Universe. The amplification {lQM) ^ Bertschinger (19g5)) Teyssier et aL (1997) and 

of primordial density fluctuations by gravity is thought to Subramanian (200 0). Cold Dark Matter (CDM) studies 

be the origin of structure formation, however the details , , , . , , ,. AT u j • w t j 

° ' based on high-resolution W-body simulations performed 

of the formation process are not yet well understood and by NavarrQj Frenk & white (1996j 1997) suggest a cuspy 

the study of the structure and properties of dark matter and universal dark matter (DM) density proffle in galax _ 

halos and of the intra cluster plasma in virialized systems ieg &nd dusters of galaxieSj independent of mass scale 

can give important clues to understand the physics in- &nd cosmology . this result is contrad i cte d by Jing & Suto 

volved m the formation and evolution of galaxy clusters. rtnnfl ^ „ . , , , , . , r , 

° J (2000). However, some important observational facts seem 
Nowadays, these studies have undergone great improve- , , , jjt,j.u j. j- • i • i 

J ' ° f * no t i 5 C reproduced by these studies: numerical Simula- 

ments with the developement of advanced observational tiong baged on thg CDM scen&rio predict density profiles 

facilities and techniques, together with the progress of nu- ^ gteep inner glopes which m tQ reproduce the rotation 

merica simu a ions. curves of low surface brightness (LSB) galaxies (Flores & 

Many works have been developed during the last Primack 1994; Moore et al 1999) _ Although other wor ks 

decades on this respect. Secondary mfall and the ef- dgim ^ cugpy DM profiles are consistent with the avail _ 

fects of this process o n the cluster structure were dis- abk dat& for dwarfg and LgB galaxieg (yan den Bogch & 

Send offprint requests to: R. Demarco (demarco@iap.frj Swaters 2001), microlensing studies towards the center of 

* Based on data obtained from the public archive of ROSAT our galaxy also support the incompatibility between CDM 

at the Max-Planck-Institut fur Extraterrestrische Physik. simulations and observational evidence (Binney & Evans 
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2001). To explain the discrepancy concerning the central 
slopes of DM halos (Navarro, Frenk & White 1997; Moore 
et al. 1999), the lack of sufficient resolution in the cen- 
tral regions of simulated halos has been proposed; Point 
Spread Function (PSF) effects together with insufficient 
resolution on galaxy rotation curves (van den Bosch & 
Swaters 2001) may be also at the origin of the discrep- 
ancy between models and observations. 

Observations in the X-ray band provide valuable in- 
formation on the hot Intra Cluster Medium (ICM). A 
popular model used to fit the spatial distribution of the 
X-ray gas is the so called /3-model (Cavaliere & Fusco- 
Femiano 1976; Sarazin 1988) which assumes an isother- 
mal ICM. However, this model may not be good enough 
to describe the cluster gas component, since the isother- 
mality of the ICM is still rather controversial. Cooling 
flows are known to produce a drop of the gas tempera- 
ture towards the center of the cluster; besides, outside the 
cooling flow region things are still not clear. Temperature 
profiles based on ASCA (White 2000) and ROSAT (Irwin, 
Bregman & Evrard 1999) data were found to be consis- 
tent with the isothermal hypothesis. On the other hand, 
Markevitch et al. (1998) found from ASCA observations 
that cluster temperatures decrease significantly with ra- 
dius. Irwin & Bregman (2000) analysed BcppoSAX data 
and claimed a slight rise in cluster temperature with ra- 
dius, a result which is at odds with De Grandi & Molendi 
(2002), who found that for a set of 21 clusters the tempera- 
ture profile has a clear isothermal core (excluding the cool- 
ing flow region) followed by a rapid radial decline. Note 
that such a core is consistent with Chandra observations 
of cooling flow clusters, where the temperature profile rises 
rapidly with radius, then remains approximately constant 
out to ~ 0.8 Mpc (Allen, Schmidt & Fabian 2001). Finally, 
XMM-Newton observations e.g. of Coma (Arnaud et al. 

2001) and Abell 1795 (Tamura et al. 2001) show small but 
significant radial variations of the temperature. Numerical 
simulations seem to confirm a decline of temperature pro- 
files with radius, but are not able to reproduce the flat- 
ness of these profiles in the innermost regions (Frenk et al. 
1999; Loken et al. 2002). This disagreement may be due 
to additional physical processes that must be taken into 
account in future numerical simulations. Another point is 
that a single (3 value cannot always fit the X-ray surface 
brightness profile of clusters (Allen, Ettori & Fabian 2001; 
Allen, Schmidt k Fabian 2001; Hicks et al. 2002). 

The ICM density and temperature distributions are 
of fundamental importance because they can be used to 
determine the specific entropy distribution of the ICM, 
thus providing important information to understand non- 
gravitational internal and external processes that may 
contribute to the ICM thermal history, such as external 
preheating and energy injection supernova-driven galaxy 
winds (Brighenti & Mathews 2001; Dos Santos & Dore 

2002) . Non-gravitational processes may be responsible for 
the observed breaking of the self-similar relation between 
X-ray luminosity and temperature predicted by theory 
(Arnaud & Evrard, 1999; Rosati, Borgani & Norman 2002 



and references therein) and also the so called Entropy 
Floor (Ponman, Cannon & Navarro, 1999; Helsdon & 
Ponman 2000; Lloyd-Davies et al. 2000). 

With the assumption that the X-ray plasma in clus- 
ters of galaxies is weakly dissipative, clusters considered 
as self-gravitating systems are likely to verify properties 
similar to those recently found in elliptical galaxies, con- 
sidered as self-gravitating systems. Namely, the optical 
surface brightness profiles of elliptical galaxies can be fit 
by a Sersic law (Sersic 1968; Caon et al. 1993; Ciotti & 
Bertin 1999): 



S(s) = Sq exp 



(1) 



characterized by three parameters: So (intensity), a (scal- 
ing) and v (shape). For a sample of 132 ellipticals be- 
longing to three galaxy clusters, the Sersic parameters 
were found to be correlated two by two, and in the three- 
dimensional space defined by these three parameters they 
are located on a thin line. These properties have been 
interpreted as due to the fact that, to a first approxima- 
tion, all these elliptical galaxies have the same specific 
entropy (entropy per unit mass) (Gerbal et al. 1997, Lima 
Neto et al. 1999, Marquez et al. 2000), and that a scaling 
law exists between the potential energy U and the mass 
M for these galaxies: U cx Af 1 - 72 ^ ' 03 (Marquez et al. 
2001). Each of these relations defines a two-manifold in the 
[ZogSo, loga, v] space. The thin line on which the galaxies 
are distributed in this space is the intersection of these 
two two-manifolds. Such relations are most probably a 
consequence of the formation and evolution processes un- 
dergone by these objects, since theory predicts U cx Af 5 / 3 
under the hypothesis that energy and mass are conserved 
(Marquez et al. 2001). 

Interestingly, numerical simulations of cold dark mat- 
ter haloes in two different mass ranges lead to a similar 
scaling law between the potential energy and mass of the 
haloes. In the mass range 4 x 10 5 < M < 4 x 10 8 M Q 
(unvirialized clusters), Jang-Condell & Hernquist (2001) 
find a relation consistent with U cx A/ 5 / 3 , while in the 
mass range 10 12 < M < 10 15 M Q (virialized clusters) 
Lanzoni (2000) finds U cx M l l39±0 02 . 

In this work, we present a study aimed at testing 
whether results similar to those found in elliptical galax- 
ies can also be obtained for galaxy clusters, based on an 
accurate modeling of the cluster X-ray surface brightness. 
We use a de-projection of the Sersic profile (Eq. to ob- 
tain the gas and DM density distributions, temperature 
profiles, dynamical mass distributions and estimations of 
the integrated specific entropies of the gas and DM com- 
ponents for a set of 24 nearby galaxy clusters and a group. 
Interesting correlations between physical quantities are 
found, comparable to those observed in elliptical galax- 
ies, which can give important clues to understand better 
the formation and evolution of galaxy clusters. This paper 
is structured as follows: our sample is described in Sect.|2 
the calculations of the physical quantities used in this pa- 
per are presented in Sect. EH the method used to determine 
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the gas density profile from the X-ray surface brightness 
is described in Sect. ^ the methods used to derive the 
temperature profile and the dark matter distribution are 
explained in Sect. El results are presented in Sect. El and 
conclusions in Sect. El 

2. The sample 

We have retrieved data taken with the PSPC-B camera 
of ROSAT from the ROSAT archive at MPE. The en- 
ergy range considered is 0.44—2 keV, corresponding to 
the four energy bands R4 to R7 (Snowden et al. 1994). 
These bands were chosen in order to avoid the low signal 
to noise ratio in the lower bands due to the high absorp- 
tion by the hydrogen column. The spatial resolution of 
the PSPC is 25" and its energy resolution corresponds 
to 0.43% at 0.93 keV. We selected observations with the 
longest exposure times and where the cluster showed a reg- 
ular shape, with no obvious mergers and a smooth light 
curve (no strong scattered solar X-ray contamination) . We 
thus built a sample of 24 clusters (see Tabled with red- 
shifts ranging between 0.01 and 0.3. Redshifts were taken 
from the SIMBAD data base (except for A2199 for which 
the redshift was obtained from Wu, Xue & Fang 1999), and 
gas temperatures and luminosities from Wu, Xue & Fang 
(1999), except for A2034 and A2382 for which tempera- 
tures were taken from Ebeling et al. (1996). We assume 
H = 50 km s" 1 Mpc" 1 , Q = 1, and A = through- 
out this analysis. In order to increase the range in Tx, 
in particular to include cooler systems when drawing the 
Lx — Tx relation, we intended to add several groups to 
our sample of clusters. However, we only included in our 
sample the group HCG 62, the one with the longest expo- 
sure time and best signal to noise ratio (we also tried to 
include HCG 94, but discarded it because of its signal to 
noise ratio). More groups will be considered in forecoming 
works. 

The data reduction was done using the software devel- 
oped by Snowden et al. (1994). The routines in the soft- 
ware provide the best available modeling and subtraction 
of various non-cosmic background components and cor- 
rections for exposure, satellite wobbling, vignetting and 
variations of detector quantum efficiency. A flat-field cor- 
rection of the images was applied and the non-extended 
sources were masked, except the cluster centers. The whole 
procedure was carried out only in clusters without strong 
scattered solar X-ray contamination; for each cluster, we 
checked the light curves in the 4 energy bands considered, 
and all those with count rate peaks larger than 3 counts 
sec -1 in their light curves were excluded. 



3. Estimating physical quantities 

3.1. Gas density profile 

The observed X-ray emission of the ICM is directly related 
to the gas distribution in the dark matter halo gravita- 



tional potential. Thus, in order to compare theory with ob- 
servations, a description of the gas distribution is needed. 
Using a parameter-dependent model for the gas density 
profile, it is possible to re-construct the 3D X-ray emis- 
sion of the cluster which, once projected and compared 
to the observations (Sect. 0J|, will allow us to derive the 
best set of values for the model parameters. We have cho- 
sen a 3D deprojection of a Sersic profile (Sersic 1968) to 
describe the gas distribution in clusters. This choice was 
motivated by the fact that we already used this profile to 
fit the optical surface brightness of elliptical galaxies, and 
computed all the physical quantities needed here, such as 
the entropy, potential energy, etc., as a function of the 
three Sersic parameters (see Marquez et al. 2001 and ref- 
erences therein). Note that from a mathematical point of 
view, since the Sersic profile has three parameters instead 
of two (compared to other models as for instance the (3- 
model), the fitting process is more flexible. Besides, the 
fact that the volume integral of this profile does not di- 
verge at large radii allows us to compute important quan- 
tities such as the total mass, potential energy and entropy 
of the system without any extra mathematical require- 
ment such as a cutoff radius, for instance. Note also that 
the Sersic law (Eq. (JTJ) is a non- homologous generalization 
of the de Vaucouleurs i? 1 ' 4 profile (de Vaucouleurs 1948). 
The 3D deprojection of such a profile corresponds to a 
generalized form of the Mellicr-Mathez profile (Mcllier & 
Mathez 1987) given by: 

Pgas{r) = po {r/ay p exp[-(r/a) v ] (2) 

where po is the volume gas density associated to the cen- 
tral column density So and the parameters p and v are 
related by the numerical approximation (Marquez et al. 
2001): 

p~ 1.0 - 0.6097^ + 0.05563i/ 2 (3) 

which gives the best approximation to the Sersic law when 
Eq. J5J) is projected. The Sersic profile defined by Eq. i ffl l 
corresponds to a surface mass density while Eq. (0 is the 
volume mass density. The condition that the mass ob- 
tained by integrating Eq. must be equal to the mass 
obtained by integrating Eq. implies: 



where T(a) is the complete gamma function defined by 
T(a) = f Q °° x a - x e- x dx. 

3.2. Dark matter distribution and dynamical mass 

Once the gas distribution is known, a reasonable hypoth- 
esis can be used to derive the dark matter distribution 
in the cluster. Previous works on X-ray clusters suggest a 
power law relation between the distributions of dark mat- 
ter and gas (e.g. Gerbal et al. 1992; Durret et al. 1994). 
We will assume here a relative distribution of the DM and 
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gas of the form Pdm/ Pgaa = R{ r )i where R(r) is a power 
law of the form: 



R(r) = K (- y 



(5) 



Under this hypothesis, the dark matter also follows 
a Sersic law: it decreases exponentially above a certain 
radius and its asymptotic behavior towards the cluster 
centre is a power law of slope p' = —(p+ot). The values 
for k, assuming a = 0.25 are given in Table[5](see Sect.0J|. 

Using Eqs. JSJ and (J2J), the total amount of mass con- 
tained within a spherical region of radius r is given by the 
integral: 



p = 0.6 for the ICM) and m p the proton mass. The elec- 
tron number density and the gas mass density are related 
by n gas (r) = p gas (r)/(1.14 m p ). The gas temperature pro- 
file can be obtained as a function of k and a by replacing 
Eqs. 10 and © into Eq. and performing a Gauss- 
Laguerre integration of the latter. The solution is of the 
form: 



T(r,K,a)—( — \ (-\ e(") t(r,K,a) 



M Dyn (r, k, a) 







h 1 


f 




/o 







k 7 



3 — (p + a) (T\ v 



■(;) 



+ 7 



p gas (u)A-KU 2 au 
3 — p fr\ u 



(6) 



where 7(0,2) is the incomplete gamma function defined 

by 7(0, z) — L x a ~ 1 e~ x dx. 

3.3. Gas temperature profile 

An important point in our study is to compute the tem- 
perature distribution of the ICM. This can be achieved 
by estimating the gas density profile, obtained by fitting 
the observations, and assuming that clusters of galaxies 
are systems in a nearly hydrostatical equilibrium sate. A 
hypothesis on the state equation of the ICM gas is also 
needed. An ideal gas state equation can be considered as 
a good approximation, although its application to self- 
gravitating systems has been questioned in the past (see 
Bonnor 1956 and references therein). 

Therefore, the equation of hydrostatical equilibrium: 



VP(r) = -G%W Ur) 



(7) 



is then combined with the equation of state for the hot 
intra-cluster plasma: 



P{r) = ^^k B T gas (r) 
pm p 



(8) 



to provide the following equation from which the ICM 
temperature as a function of radius T gas (r) can be derived 
once the gas number density as a function of radius n ga s(^) 
is known: 



+ 



dT gas (r) 
dr 



(9) 



where G is the gravitational constant, ks is the Boltzman 
constant, p is the plasma molecular weight (we assume 



(10) 



where w = 4ttG ^ ^- = 1.54 x 10 38 ^ m keV, 
and 

3 — (p + a) 



(11) 



t(r,K,a) = I <; K7 



3-p 



-,x 



(p+*+i) __ , 
x e ax 



Eq. IjlOl) will be used together with a model for the hot 
plasma spectrum to derive the best set of values for n and 
a, in agreement with quantities already observed in clus- 
ters of galaxies such as their mean X-ray temperature and 
luminosity and their baryonic mass fraction (see Sect. [5j) . 

3.4. Potential energy and specific entropy 

Changing the upper limit of the integral in Eq. Ip?)l to 00 
we obtain the total dynamical mass: 



M. 



Dyn 



Anpoa 3 



3 - (p + a) 



r 



3-P 



(12) 



The total potential energy of the cluster is given by: 
MP 



Upot — G- 



Dyn 
Ra 



(13) 



where the gravitational radius R g is defined by R g = aR*, 
where a is the scale parameter and R* is a dimensionless 
radius given by the numerical approximation: 



ln(R* Q 



0.82032 - 0.92446 ln(i/) 



0.84543 



(Marquez et al. 2001). 

In spite of the X-ray emission, which is responsible in 
many cases for cooling flow processes affecting the equi- 
librium state of the cluster in the inner regions, we may 
consider clusters as structures where dissipating processes 
are negligible compared to their gravitational energy, thus 
settling into a nearly thermodynamic equilibrium at large 
scale. This can be inferred by a simple order of magni- 
tude calculation: the potential energy of a cluster of mass 
~ 10 15 Af Q and radius - 1 Mpc is about 8 x 10 64 ergs. 
The energy lost through X-ray emission during a Hubble 
time (~ 2x 10 17 s) by a cluster of X-ray luminosity ~ 10 45 
erg s _1 is around 3 x 10 62 ergs, a value which is almost 300 
times smaller than its potential energy. We can therefore 
estimate the gas entropy of such a configuration. 
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The specific entropy of the intra-cluster gas can be 
computed from the distribution function in the phase 
space, /(x, v), of the gas particles using the microscopic 
Boltzmann-Gibbs definition: 



Jfln(f) d 3 xd 3 v 
J f cPx <Pv 



(14) 



where the Boltzmann constant is ks = 1. Note that this 
expresion gives us the specific entropy of the entire system 
(gas or DM) because the integration covers the total vol- 
ume in phase space. This definition then corresponds to 
an integrated specific entropy, and it will be referred to as 
'global' specific entropy. It is important to say that when 
we use the words 'integrated' or 'global' for the specific 
entropy, we are referring to this definition applied to the 
gas or DM separately and not to the sum of these two 
components. 

To find the distribution function some simplifying hy- 
potheses are needed. The first one is that our system is 
spherically symmetric, and the second one is that the ve- 
locity dispersion of the gas particles is isotropic (we ne- 
glect any possible rotation of the gas). Thus the distri- 
bution function /, depending explicitly only on the total 
energy, can be obtained from the density profile by an 
Abel inversion (Binney & Tremaine 1987). In this way, 
/ = /(/?) and the integrated specific entropy can be com- 
puted numerically as a function of the Sersic parameters 
only, providing a way to estimate the integrated specific 
entropy of the ICM from its surface brightness fit. It is 
important to emphasize here that the Boltzmann-Gibbs 
definition of the specific entropy makes no assumption on 
the equation of state, in particular the ideal gas equation 
for the ICM, and no assumption either on the validity of 
a hydrostatic equilibrium state of the system; this defini- 
tion is therefore a general one since it only assumes spher- 
ical symmetry and an isotropic velocity dispersion tensor. 
Since the method we present here, based on the cluster X- 
ray surface brightness fitting, provides a good constrain 
on the gas density profile, this density profile, together 
with Eq. H14|) can be used to obtain a good estimate of 
the global ICM specific entropy. 

The specific entropy can also be obtained from the 
following set of equations (Balogh, Babul & Patton 1999): 



s(r) 



and 



K {r) 



2v(pm p f/ 3 K {r) 
h 2 



kBTga S {r) 
pm p p gas (r) 2 / 3 



(15) 



(16) 



where c v is the specific heat capacity at constant volume 
of the plasma. These widely used equations require, how- 
ever, the following assumptions on the ICM : the gas is 
considered to be monoatomic and the equation of state 
for an ideal gas is supposed to hold. By knowing the ICM 
density and temperature distributions, Eqs. I|15|) and (|16|l 



can then be used (see Sect. 16. 311 to compute the gas spe- 
cific entropy profile for each cluster. Gas density profiles 
are obtained from the X-ray surface brightness fit (see 
Sect. while temperature profiles are derived by assum- 
ing both the hydrostatic equilibrium condition and the 
ideal gas state equation for the ICM, as explained before. 
Since the physical quantities involved here are intensive, 
these equations cannot be integrated to derive the inte- 
grated specific entropy of the gas, and Eq. (|14f> will have 
to be used for this purpose. 

The entropy is of fundamental importance to under- 
stand the effects of non-gravitational processes on the 
thermodynamical history of clusters. Previous studies re- 
fer to the gas entropy at 0.1r„j r , where r V { r is the clus- 
ter virial radius (e.g. Ponman, Cannon & Navarro, 1999; 
Lloyd-Davies, Ponman & Cannon 2000), while in this work 
we estimate the gas specific entropy for the entire clus- 
ter. Our calculations take into account all possible sources 
of heating and cooling, regardless of the distance to the 
centre, thus providing a good quantitative basis to which 
models can be compared in order to disentangle the dif- 
ferent processes that affect the internal energy of the ICM 
during its history since the earlier epochs of the cluster 
formation. 

We also notice that our assumption concerning the 
Pdm/ ' Pgas ratio (Eq. JSJ) implies that the resulting DM 
density distribution is also Sersic-like, with a central den- 
sity distribution described by a power law varying as 
r -(p+a)_ Tj n d er the hypothesis already stated concerning 
the distribution function of particles in phase space, we 
can in principle also compute the integrated specific en- 
tropy of the DM halo by using Eq. (fT4"|) and the DM den- 
sity profile. Therefore the DM specific entropy is calcu- 
lated numerically from the equation: 

1 /"* (0) dM 

' -lnf(e)de, (17) 



M 



Ac 



which can be derived from Eq. i|14|) , and where e represents 
the binding energy, M the total halo mass, and '3/ the 
relative gravitational potential (see Magnard 2002). 

4. Fitting method and density profile 

Finding the gas density profile amounts to deriving the 
best set of values for the Sersic parameters. This has 
required to fit the ROSAT images by a pixel-to-pixel 
method, which creates a three-dimensional model of the 
X-ray emission which is then projected by integration 
along the line of sight, taking into account the energy re- 
sponse and the point spread function of the detector. The 
result is a synthetical image which can be compared to the 
observation, and the best set of Sersic parameters is ob- 
tained by successive iterations. The gas density profile de- 
fined by Eq. (J2J) is used to model the bremsstrahlung emis- 
sion, and the free-bound and bound-bound X-ray emis- 
sions are taken into account. The code computes the X- 
ray emissivity e v in every point; it is then projected to ob- 
tain the surface brightness. This projected image is then 
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convolved with the ROSAT point spread function (PSF), 
which varies as a function of position (and energy) on the 
detector. We have used in our fits a FWHM of 2 pixels 
which corresponds to the central PSPC PSF. 

The cluster redshift and the mean gas temperature are 
required as input parameters. To obtain an initial guess for 
the free parameters in Eq. Q we fit a Sersic profile (Eq. 
|QJ) to the X-ray surface brightness of each cluster. In or- 
der to stay as close as possible to our hypothesis of spher- 
ical symmetry, we selected clusters presenting the most 
round and uniform projected X-ray emission. However, 
this emission is not perfectly circular and due also to the 
fact that what we observe is always a projection on the 
plane of the sky of a three dimensional structure, we de- 
cided to take into account during the fitting process the 
cllipticity of the projected emission supposing also that 
all the clusters are oblate spheroids, i.e., that the third 
axis along the line of sight is a major axis. If a is the 
semi-major axis, which also corresponds to the scale pa- 
rameter in the Sersic profile, and b is the semi-minor axis, 
we define the ellipticity of the projected X-ray emission 
by e = yJl — (b/a) 2 . The ellipticity and semi-major axis 
position angle of the X-ray distribution are thus consid- 
ered as free parameters for the fit and given as inputs for 
the code. The gas was assumed to be isothermal as a first 
approximation. 




23 h 26 m 20 s 26 m e 40 s 20 s 25 m 3 40 s 20 s 



Fig. 1. Comparison of our best model of the X-ray sur- 
face brightness (dot-dashed ellipses) with the correspond- 
ing ROSAT image for the cluster A2597. The X-ray iso- 
contours are 0.5, 1, 3, 10 and 100 times the background; 
they were made after smoothing the original image with 
a gaussian kernel of a — 3.2 pixels. 



After obtaining the initial guesses for each free param- 
eter, we use these values together with Eq. to make 
a new synthetic image, then compare it to the actual 
ROSAT image. The parameter values are then changed 
and the comparison process is repeated iteratively, until 
it finds the 3D X-ray emission which best fits the surface 
brightness profile of the observation when projected. The 
fitting process is carried out with the MIGRAD method in 
the CERN MINUIT library (James 1994). In this process, 
the gas temperature is kept constant, as a first approxi- 
mation. The exact shape of the temperature profile is not 
crucial when finding the set of values for the parameters, 
because the emissivity does not have a strong dependence 
on temperature as compared to its dependence on den- 
sity; indeed, we find that modifying the temperature pro- 
file changes the values of the Sersic parameters by at most 
a few percent, as explained in Sect.[5J 

The fitting process directly gives us the best estimates 
for the semimajor axis a and shape parameter v. However, 
the central electronic number density n e0 given by the fit is 
not accurate enough, and is estimated by another method 
(see Sect. . 

We evaluated errors with the MINUIT error function 
MINOS, which calculates parameter errors taking into 
account both parameter correlations and non-linearities. 
Resulting errors correspond to a la deviation. 

All the equations presented in Sect. |3| refer to spher- 
ical symmetry, for which we had to define an equivalent 
scale parameter, a eq in order to go from the oblate geom- 
etry considered in the fit to the spherical geometry of the 
model. This new scale is defined as a ea = (a 2 b)^ 3 and will 
be used instead of a to compute the specific entropy, the 
dynamical mass and the potential energy of a spherically 
symmetric X-ray region. 

We show in Fig. ^ a comparison of our best fit Sersic 
model of A2597. The good fit of the cluster surface bright- 
ness can clearly be seen, confirming the capability of the 
Sersic profile to reproduce the cluster X-ray emission. 



5. Gas temperature and dark matter to gas ratio 

The fitting process described above provides the best 
values of a eq and v for a given gas temperature pro- 
file, assuming as a first approximation that each clus- 
ter has an isothermal ICM, and then obtaining, through 
an iterative procedure, the best compatible gas density 
and non-isothermal temperature profiles, assuming hydro- 
static equilibrium and spherical symmetry. 

In order for the n e0 normalization to give the observed 
number of counts, we first computed the cluster flux F, 
after substraction of point sources and background (ob- 
tained by the 2D fit). The masked source pixel counts 
were set to the mean value taken within ellipses. Then 
n e0 was chosen so that a bremsstrahlung XSPEC model 
convolved by a hydrogen column absorption would give 
the same flux as seen by ROSAT: 
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F = 



3.02 x 10 



-15 



n e njdV 



3.02 x 10~ 15 2 2+2E Tr L 
~^d\ 



3-2p 



where is the luminosity distance to the cluster, and n/ 
is the number density of plasma ions. 

Eqs. Ill0[l and (f lift provide the ICM temperature pro- 
files as families of solutions depending on the parameters 
k and a. Since these are found to be different from those 
assumed to estimate the Sersic parameters, the fitting pro- 
cess must be repeated, using the new temperature profile, 
in order to obtain a new set of parameters for each cluster, 
until convergence. 

We tested on one cluster (A2029) the effect of a tem- 
perature gradient on the gas density profile from distances 
close to r e ff/4 outwards, where r e // corresponds to an ef- 
fective radius which contains half of the cluster gas mass 
(see Table EJ. For this, we did a second fit of the sur- 
face brightness of A2029 using a new temperature profile 
given by a power law of the form T(r) = To(r /r e f f)@ 
where r e ff ~ 2800 kpc. To was set in order to have a 
non-weighted mean temperature equal to the mean clus- 
ter temperature (Table ^) and we considered the cases 
j3 = —0.5 and —1. After the first iteration, the values 
found for n e0 , a eq , and v remained almost unchanged for 
both values of f3: the scale and shape parameters changed 
by about 1%, and the central electronic density by 4% 
with respect to the isothermal fit. 

We therefore decided to keep the Sersic values given 
by the original fit as the definitive ones. Eq. (|10|l can be 
re- written in the form T(r, k, a) — K,T\{r, a) + l2(r). We 
produced a set of Ti (r, a) profiles corresponding to val- 
ues of a between and 2 with steps of Aa = 0.1. These 
curves together with their corresponding T2 (r) profiles and 
an adequate set of values for k were used to find which 
combination of k and a gives the closest values to the ob- 
served global temperatures and luminosities (Ebeling et 
al. 1996; Wu, Xue & Fang 1999), and baryon mass frac- 
tions (Arnaud & Evrard 1999; Mohr, Mathiesen & Evrard 
1999; Schindler 1999). These quantities are evaluated as 
follows. The two profiles n e {r) and T(r, k, a) plus the hy- 
pothesis of a gas metallicity equal to 0.3 Z@ (see Renzini 
1997) and the emission process determine the cluster X- 
ray emission. We used the XSPEC software to simulate 
the corresponding spectra with a bremsstrahlung emission 
model (to be coherent with the emission model used in the 
fitting program) . We summed up the spectra from shells of 
constant electronic density and temperature over radius, 
considering only the volume intersected by the ROSAT 
observation cylinder. The spectrum obtained is then con- 
volved with a photoelectric absorption model and with the 
ROSAT response function to produce the simulated spec- 
trum. A fit is then performed on this spectrum to derive 
the X-ray temperature and luminosity (the metallicity and 
hydrogen column density are fixed). 



The set of n and a has to produce Lx and T in agree- 
ment with the observed data within errors. Moreover, the 
gas mass fractions should stay inside the observed limits. 
So we performed a minimization of the distance between 
the observed data and the predictions from our model. We 
found that the best a is often very close to zero, which is 
not in good agreement with results obtained from numer- 
ical simulations by Teyssier (2002) which show a pdm/ Pg 
ratio varying globally as r -°- 25 ; besides our values of a 
are not well constrained. We therefore chose to impose 
a = —0.25; n was then recomputed to give luminosities 
and temperatures as close as possible to the observations. 

To each temperature profile T(r, k, a) corresponds a 
couple of simulated observational parameters (T s i m , L s im). 
The value of n is constrained by imposing these parame- 
ters to be close (within error bars) to the real observational 
values. The gas mass fraction is checked to be compatible 
with the isothermal hypothesis. 

6. Results 

6.1. Mass distributions and parameter correlations 

The 3D gas density profiles were computed with Eq. J2Jl 
for all the clusters in our sample, using the sets of param- 
eters obtained from the surface brightness fitting process 
(described in Sect. 0} of the PSPC images and listed in 
Table|21 By means of Eq. JSJ , the corresponding DM distri- 
bution can be recovered. In Fig. [2] we show the 3D gas and 
DM density profiles for every cluster divided by the critical 
density of the Universe at the cluster redshift (e.g., Ettori, 
De Grandi & Molendi 2002): p c (z) = 3H 2 {z)/(8ttG), 
where the Hubble constant at redshift z is defined by 

h(z) = H ^n m (\ + zf + i-n m . 



Both DM and gas distributions are Sersic-like, imply- 
ing that the corresponding profiles decrease exponentially 
outwards above a certain distance from the cluster cen- 
tre; on the other hand, when r goes to zero they follow 
a power law in radius with a logarithmic slope tending 
asymptotically to —p{v) for the gas and to —{jp(v) + a) 
for the DM. 

The values for k, obtained as described in Sect. |21 are 
in most cases smaller than 10, while the value for a was 
fixed to 0.25 (see above). Although DM halos are denser 
than the ICM gas, as implied by the k factor and the 
asymptotic slope difference a, the general shapes of both 
profiles look quite similar, implying that dark matter and 
gas are distributed in a comparable way. This effect is a 
natural consequence of the hypotheses and conditions im- 
posed to our model. However such a behaviour seems to 
be the case for massive systems (like those in our sample - 
see total masses in Table 01, as shown by observations of 
galaxy clusters at moderate and high redshifts (Schindler 
1999). The gas would be accreted into the forming struc- 
ture and once the system reaches a relaxed state, the gas 
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Fig. 2. 3D deprojected gas (full lines) and dark matter 
(dashed lines) density profiles divided by the critical den- 
sity of the Universe at the cluster redshift. 



just accommodates into the halo potential. This is true 
at the scale of massive galaxy clusters where the poten- 
tial well is deep enough to prevent the gas from expansion 
due to non-gravothermal processes. In this case, both dark 
matter and gas present similar distributions in contrast 
with what is observed in smaller systems such as groups 
of galaxies and even galaxies, where the gas can produce 
extended cores as result of energy injection due to super- 
nova explosions or shock winds (Bryan 2000, Bower et al. 
2001, Brighenti & Mathews 2001). It is also important to 
mention that Eq. ©, based on galaxy cluster observa- 
tions, may no longer be valid for low mass systems such 
as groups of galaxies, in which case our DM model would 
be inappropriate to describe groups. 

The averaged gas density profile for our set of 24 galaxy 
clusters is well fit by a Sersic profile with parameters: 
po = 7.4 x 10~ 24 kg m~ 3 , a eq = 367 kpc and v = 0.56. The 
latter gives p{v) = 0.67 which makes the corresponding 
DM density profile vary as r -0 92 when r goes to zero. This 
central slope for the DM halo is shallower than the self- 
similar solution for spherical collapse in an expanding uni- 
verse found by Bertschingcr (1985) (p oc r - 2 - 25 ) ; than the 
asymptotic behaviour found by Moore et al. (1999) in their 
numerical simulations (p oc r -15 ) and than the NFW 
(Navarro, Frenk & White, 1996, 1997) universal density 
profile. However it is steeper than the p oc r~ 75 critical 
solution found by Taylor & Navarro (2001) for galaxy- 
sized CDM halos based on the study of their phase-space 
density distribution. This critical solution can be inter- 
preted as a maximally "mixed" configuration, where the 
phase-space distribution across the system is the most uni- 
form one compatible with a monotonically decreasing den- 
sity profile and with the power-law entropy distribution. 



This configuration would be the result of non-spherically 
symmetric formation processes through hierarchical merg- 
ing. Mass shells are continously mixed and the density 
profiles tend to be shallower than the NFW profile at the 
center, converging to the p oc r -0 75 distribution for max- 
imal mixing (Taylor & Navarro 2001). 

One important advantage of the Sersic density profile 
is that its volume integral converges when integrated up 
to infinity, making it thus possible to determine the total 
cluster mass, and the potential energy and entropy of the 
system without introducing a cutoff radius, in contrast to 
the popular /?-model, for instance. Total dynamical masses 
can thus be computed (see Eq. (|12|0 and the resulting 
values for all clusters are given in Table EH 




log(r/r aoo ) 

Fig. 3. Synthetic cumulative dynamical mass profiles nor- 
malized to the corresponding total integrated cluster mass 
M-Byn- The radial coordinate is normalized to the cluster 
r 2 oo radius. 



Fig. El shows the resulting cumulative dynamical mass 
profiles for every cluster as obtained by means of Eq. 10 
with the corresponding Sersic parameters. Every profile 
has been normalized to the corresponding cluster total 
mass, Mdyn (see Eq. I|12|) 'l. and the radial coordinate has 
been normalized to the corresponding cluster r2oo radius, 
which is defined as the radius within which the mean den- 
sity is 200 times the critical density of the universe, p c , as 
defined above. In general, we can define a radius ta within 
which the mean density is A times p c and ta can be ob- 
tained from the relation M(r < rA)/(47rr A /3) = Ap c ). 
The averaged cumulative dynamical mass distributions 
show a logarithmic slope dlogM / ' dlogr ~ 1 at r ~ 0.7r2oo- 
Further out this slope decreases rather fast. According 
to our model, the cumulative mass profiles converge only 
around r ~ 1 Or 200, but the mass variation is only of a 
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factor of two, going from M ~ 10 15 Mq at ~ r2oo to 
M~2x 10 15 M Q at - 10r 200 - 

The analysis by Ettori, De Grandi & Molendi (2002) 
of BeppoSAX data includes some clusters of our sample. 
They estimate masses for A = 1000 and 2500, by using 
either a NFW or a King profile for the total mass distri- 
bution. Comparing our results with theirs for same clus- 
ters, we see that our mass estimates, derived from the 
Sersic profile for these two values of A within same ta 
radii as indicated in their Table 2 are in quite good agree- 
ment. In most cases, our masses differ by about 10% or 
less, with differences reaching about 20% for a few cases. 
These differences are likely to be due to the use of differ- 
ent functional forms (the Sersic, NFW and King models) 
to describe the mass profile. Moreover, from our sample 
we find in average r2oo ~ 2.5 ± 0.4 Mpc (see Table |2J). 
These values are about 78% of those found in numerical 
simulations by Navarro, Frenk & White 1996 for cluster 
sized systems of comparable masses. Moreover, we note 
that the values of our scale parameter a eq are compara- 
ble to those of the scale radius r s of the NFW profile 
for similar mass ranges. In this way, we obtain in aver- 
age a eq /r2oo ~ 0.12 in good agreement with the value of 
r s /rwo ~ 0.14 (Navarro, Frenk & White 1996), the dif- 
ference being of only 14%. We can therefore say that the 
Sersic profile is able to describe the mass distribution in 
clusters in as much the same way as other models as the 
NFW and King models, the differences being due to the 
mathematical natures of the models. 

Based on a de-projected Sersic model for the gas den- 
sity profile, the best set of values of n e0l a eq and v for 
each cluster is given in Table El These three parameters 
are displayed two by two in Figs.0J|Hland|ni A clear corre- 
lation between the shape (y) and scale (a eq ) parameters is 
seen, while the correlations we find for So — v and So — a eq 
(So being obtained from po by using Eq. are clear but 
show a somewhat larger scatter. Note that these three cor- 
relations have shapes similar to those found in elliptical 
galaxies. This may indicate, as for ellipticals, the existence 
of an entropic line on which galaxy clusters lie, in which 
case the correlation shown in Fig. 0] may just be the pro- 
jection of this isentropic relation on the log(v) — log(a eq ) 
plane. In the case of elliptical galaxies, this entropic line 
was interpreted as the intersection of two surfaces in the 
[log S , log a, log v] space : the Entropic surface and the 
Energy-Mass surface (Marquez et al. 2001). In this work, 
we find that these two surfaces exist for our galaxy clus- 
ters as well, and these clusters are also located on a line 
corresponding to the intersections of the two surfaces (see 
Sect. 16.51) . A complete discussion on this point will be pre- 
sented in a forecoming paper (Magnard, in preparation). 
It is worth noting that, although from the mathematical 
point of view no correlation between the model parame- 
ters is expected, the fact that we observe such correlations 
probably indicates that they are due to the physics under- 
lying the X-ray emission distribution. 



Fig. 4. Correlation between the density profile parameters 
a ea and v. 
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Fig. 5. Correlation between the density profile parameters 
Sq and v. 
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Fig. 6. Correlation between the density profile parameters 
Sq and a eq . 
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6.2. Temperature profiles 

We show in Fig.0the 3D non- weighted temperature pro- 
file for every cluster in our sample, calculated with the 
parameters given in Table [21 and normalized to the cor- 
responding median temperature from the literature (see 
Table UJ. As for the total mass profiles, the radial coor- 
dinate has been normalized to the corresponding cluster 
virial radius. All the profiles are very similar: they slightly 
increase from the center, then remain approximately con- 
stant to finally decrease at large radii. The mean tem- 
perature profile is consistent with an almost isothermal 
distribution for radii smaller than 0.25r2oo (~ 625 kpc), 
in agreement with other works (Allen, Schmidt & Fabian 
2001; De Grandi & Molendi 2002), followed by a rapid 
decrease for radii larger than 0.25r2oo, with a logarithmic 
slope of the order of —0.6 for radii around 0.63 r2oo- This 
decrease results to be quite similar to the slope of —0.64 
found by De Grandi & Molendi (2002) at radii larger than 
~ 750 kpc for their mean temperature profile, taking into 
account the cooling and non-cooling flow clusters of their 
sample. The results of Markevitch et al. (1998) based on 
ASCA observations show similar trends for the tempera- 
ture distribution. The characteristic central drop we ob- 
tain in our temperature profiles is due to the mathematical 
properties of the density profile rather than to a real phys- 
ical effect, but the lack of resolution of the PSPC impedes 
from properly addressing this issue, and we cannot tell 
anything about the temperature distribution within the 
central cooling flow region. 



In Fig. IS] we show the re-projected emission weighted 
temperature profile for every cluster, defined by : 

J n 2 c {r)T{r)dz 

T CW {R) — — t. 2? (18) 

The Lx — Tx relation derived from our calculations is 
shown in Fig. |§| Our power law fit (excluding HCG 62) 
is Lx oc j^.- 65±0 - 17 j the exponent is in agreement with 
Markevitch (1998), who finds L x oc T| 64±0 - 27 , but our 
line is shifted towards higher luminosities. This shift is 
probably due to the fact that we include the central region 
in our luminosity evaluation, while Markevitch excises it. 

6.3. ICM and DM specific entropies 

The global specific entropy s gas of the ICM is given in 
Table It is found to vary very little from one cluster 
to another, as for the specific entropy of stars in ellipti- 
cal galaxies (Marquez et al. 2001 and references therein). 




-1.5 -1.0 -0.5 0.0 

l°g(r/r soo ) 

Fig. 7. Synthetic 3 D temperature profiles calculated from 
Eqs. i|10|) and l|ll|) with the values of a and k given in 
Table|2l The temperature is normalized to the global clus- 
ter temperature and the radial coordinate is normalized 
to the cluster r2oo radius. 
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Fig. 9. Relation between Lx and Tx- The point at the 
lower left corner corresponds to the group HCG 62. The 
solid line shows the best fit (see text) and the dotted line 
is the fit found by Markevitch (1998). 



This is, however, a first order behavior. Numerical simu- 
lations of elliptical galaxies formed in a hierarchical merg- 
ing scheme, show that their specific entropy varies a little 
with mass, most probably due to merging processes (Lima 
Neto et al. 1999, Marquez et al. 2000). The situation is not 
different in galaxy clusters. We show in Fig. ^3 that the 
global specific entropy of the ICM and the gas mass are 
indeed clearly correlated, although this is a second order 
effect (about 2%), small compared to the dominant rela- 
tion Sgas ~ constant that we observe. The difference with 
elliptical galaxies is that the slope in Fig. is steeper 
than for ellipticals: if we write s gas — sq + 0ln(M gas ), 
our best fit to the data gives — 1.86 ± 0.15 compared 
to (3 ~ 1 for ellipticals (with so constant). The error bars 
correspond to la deviations and were computed from the 
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Fig. 8. Re-projected emission weighted temperature deduced from our model. The projected radii s are normalized to 
the corresponding r v — r*2oo radii, and the temperatures are normalized to the observed isothermal temperatures. 
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Fig. 10. Relation between the gas integrated specific en- 
tropy, Sgas and the gas mass. The best fit line is indicated 
(see text). 
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Fig. 11. Relation between the gas integrated specific en- 
tropy and the mean gas temperature taken from the liter- 
ature (see text). The point at the lower left corner corre- 
sponds to the group HCG 62. 



parameter errors given by MINOS and by means of Monte 
Carlo simulations. Every parameter was modeled by a 
gaussian distribution with a equal to the corresponding 
ltr deviation from MINOS. These distributions were then 
used to derive the total mass and specific entropy distri- 
butions and the corresponding ler errors for each cluster. 
For the linear fit we used a linear least-squares approx- 
imation in one dimension, taking into account the error 
bars in both directions at the same time. A comparable 
relation is found between the integrated specific entropy 
s gas of the ICM and the dynamical mass of the cluster, 
with a somewhat steeper slope of 2.67 ± 0.12. In our fit- 
ting of the cluster X-ray surface brightness we have also 
taken into account the cluster center, thus our estimation 
of the gas specific entropy necessarily takes into account 
the effects of the central cooling flow. This information is 



integrated together with the other non-gravitational ther- 
mal processes affecting the intra cluster gas. 

The integrated specific entropy of the X-ray gas is dis- 
played in Fig. ^] as a function of the observed isothermal 
gas temperature. The entropy appears to increase with 
the temperature, consistently with a linear increase (the 
best fit to our data is shown in Fig. ^] and corresponds 
to Sg as oc T 4 - 92 ). However, the dispersion is too large to 
assert a linear dependency. Note that such a relation has 
already been observed by Ponman et al. (1999) and Lloyd- 
Davies et al. (2000), and predicted by gravitational simu- 
lations (Borgani et al. 2001, Muanwong et al. 2002). The 
only group shown in this plot does not exhibit a significant 
entropy floor, but other groups need to be added. It would 
be important to get a good evaluation of this entropy floor 
as it is a strong constraint on the non-gravitational energy 
injection (Lloyd-Davies et al. 2000). 
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Fig. 12. Kq profiles based on density and temperature 
models. 



The derived gas density and temperature profiles can 
be used to compute the gas specific entropy profile for 
each cluster by using Eqs. (|15|) and l|ltj[) . 

The Ko(r) profiles are shown in Fig.^J where the ra- 
dial coordinate is normalized to the r2oo radius. All these 
profiles are consistent with a gas specific entropy increas- 
ing with radius, indicating that the gas has been heated by 
shock fronts at successive epochs when collapsing into the 
cluster potential well, before being virialized. Thus, the 
specific entropy profile provides useful information about 
the ICM thermal history. Some of the profiles present a 
clear steepening with radius at about 0.8 r2oo, while oth- 
ers correspond to power laws throughout. The average of 
the profiles defined by Eq. I|ltj|) na s a logarithmic slope 
dlogKo/dlogr ~ 1 at about 0.6 r2oo and it is in agree- 
ment with results obtained by Tozzi & Norman (2001) 
who studied the effect of an entropy excess in the ICM 
gas before accretion into the DM halo. No entropy core is 
present, perhaps due to the effect of central cooling, al- 
though this could also be due to the mathematical nature 
of the profiles derived from the Sersic model. In general 
our profiles are in agreement with external heating models 
for rich clusters (Tozzi, Scharf & Norman 2000; Tozzi & 
Norman 2001). However, we notice that the Ko(r) profiles 
obtained from our density and temperature distributions 
are about a factor of 7 larger than those inferred by Tozzi 
& Norman (2001). This is probably just due to a normal- 
ization effect related to the way the density and temper- 
ature are computed; it is mainly the entropy variations 
that are a reliable physical quantity. 

The assumed relation ((SJ between the dark matter and 
gas densities allows us to also compute the global dark 
matter entropy (see Table EJ numerically from Eq. fTTjl . 
which comes out directly from the definition given in 




14.8 15.0 15.2 15.4 15.6 
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Fig. 14. Relation between the dynamical mass and the 
potential energy. The solid line corresponds to the best fit 
to the data and the dashed line is the theoretical slope of 
/ = 5/3 defined by the relation ln(U pot ) — I ln{Md yn ) — 
const. 



Eq. I|14[l . The corresponding entropy-mass relation thus 
found for the dark matter is very close to that obtained 
for the gas, with [3 ~ 1.90 ± 0.17 as seen in Fig. H31 



6.4. Potential energy - mass relation 

The relation between the dynamical mass and potential 
energy, previously shown for elliptical galaxies (Marquez 
et al. 2001) and in numerical simulations (Lanzoni, 2000; 
Jang-Condcll & Hcrnquist, 2001), is displayed in Fig. 1141 
for our set of clusters. If we write it as ln(U po t) — 
I In(Mdyn) = const, the best fit to our data imply / = 
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1.68 ± 0.08, in excellent agreement with the theoretical 
value of 5/3. The error bars and linear fit are computed 
in the same way as explained in the previous section. 

6.5. Discussion 

The correlations presented in the previous sections be- 
tween the mass, potential energy and integrated specific 
entropy, confirm the existence of an entropic surface, 
s(S ,a, v) — (3 ln(M(E ,a,is)) = const (see Eq. C3J), 
and a potential energy-mass surface, Zn([/ po t(So, a, v)) — 
I ln(MdynC^Oi a, v)) = const, in the Sersic parameters 
space, remarkably similar to the case of elliptical galaxies 
(Marquez et al. 2001). Indeed, we observe the clusters in 
our sample to be located at the intersection of these two 
surfaces, indicating the existence of an "entropic line" for 
galaxy clusters (see Magnard 2002). We have checked by 
simulations that this line is also recovered for a random set 
of values of the Sersic parameters, suggesting the possibil- 
ity that this relation could be a consequence of the model 
assumed to describe the final viralized system (in our case 
a Sersic profile). However, the fact that a Sersic profile 
reproduces very well the X-ray surface brightness of clus- 
ters, and hence the gas distribution of the ICM, supports 
the idea that the physical processes operating during the 
formation and evolution of galaxy clusters, which are of 
course responsible for the final structure reached by the 
ICM and DM halo, are indeed at the origin of the entropic 
line. To confirm this, the same method as in this work, but 
with different models for the gas density profile should be 
used. Moreover, the existence of both an entropic surface 
and a potential energy-mass surface for galaxy clusters 
implies that these objects can be considered as a single- 
parameter family, described by one of the Sersic parame- 
ters only (e.g. Marquez et al. 2001). Interestingly, analo- 
gous correlations have been obtained for galaxy clusters by 
Fujita & Takahara (1999). The importance of the above 
mentioned correlations resides in the fact that they are 
probably the result of the physics ruling cluster formation. 
A correlation between the global specific entropy and the 
mass conserves information on the various events affect- 
ing the thermodynamical history of clusters. The observed 
variation of s gas with dynamical mass in clusters suggests 
that dissipating processes in clusters play an important 
role as generators of entropy. These mainly correspond to 
Brcmsstrahlung emission (L oc M 4 / 3 ) and cooling flows. 
Merging processes between clusters are of importance in 
such a relation, because of their impact on the final total 
mass and on the amount of entropy produced during the 
cluster formation. Violent merger events can be accompa- 
nied by an important dissipation of energy and creation 
of entropy, while minor mergers can be translated in an 
adiabatic accretion of matter without a significant pro- 
duction of entropy. These energy losses, however, are all 
negligible compared to the cluster gravitational energy. 
Thus the value of the slope f3 in the specific entropy-mass 
relation reflects the impact of such processes on the clus- 



ter history. Clusters with higher global specific entropy 
could have undergone more episodes of hierarchical merg- 
ing through their histories, thus becoming more massive. 

On the other hand, considering the collapse of matter 
to form a virialized gravitational system, the correlation 
between the potential energy and the total mass of the 
final structure is a natural consequence of the conservation 
of energy and mass during its formation. A self-similar 
relation defined by U oc Af 5 / 3 is expected from theory 
(see Marquez et al. 2001) and we show that it is indeed 
also followed by our observed galaxy clusters. 

All these results strongly suggest that the formation 
processes affecting galaxies and clusters of galaxies are 
quite similar regardless the scale involved. 

7. Conclusions 

We have shown in the present work that the Sersic profile 
can be used as a good tracer of the matter distribution 
in clusters, under the assumption that clusters are well 
relaxed structures, as it was the case in elliptical galaxies. 
Its mathematical properties make it an interesting and 
useful tool that can be employed to explore the physics 
of relaxed systems, although any other appropriate profile 
can be used. The density profiles obtained here reproduce 
well the X-ray surface brightness profiles of the ROSAT 
PSPC images. The asymptotic behaviour of these pro- 
files towards the cluster center turns out to be shallower 
than the NFW profile, but still within the limits predicted 
by numerical simulations concerning the central slope of 
galaxy-sized DM halos. Temperature profiles derived here 
(considering the hydrostatical equilibrium hypothesis for 
the cluster structure) are in agreement with other works. 
We estimate the integrated specific entropy content for 
galaxy clusters and our specific entropy profiles are consis- 
tent with the predicted shape of the entropy distribution 
for massive clusters, obtained by simulations which take 
into account pre-heating and cooling processes. 

We have shown that both for the gas and for the dark 
matter the integrated specific entropy and the potential 
energy have behaviors comparable to those observed for 
stars in elliptical galaxies: the integrated specific entropy is 
constant to first order and in reality increases slightly with 
mass as a logarithmic function, while the potential energy 
scales with mass as a power law. Note that the index of 
this power law is close to the theoretical value of 5/3 for el- 
liptical galaxies and clusters. This strongly suggests that 
all these self gravitating systems behave similarly, even 
though they may have very different masses and thermo- 
dynamical histories. Elliptical galaxies could then be con- 
sidered as scaled down versions of galaxy clusters (Moore 
et al. 1999). Moreover, integrated specific entropy-mass 
and potential energy-mass correlations should be the re- 
sult of the formation history of the clusters. Heating mech- 
anisms and merger events play an important role here and 
total mass and energy are conserved during the whole for- 
mation process of the final virialized structure. 
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It would be interesting to apply the Sersic model to 
high rcdshift galaxy clusters in order to test the possible 
evolution of the scaling relations found in this work. The 
use of Chandra and XMM-Newton data will be crucial in 
these kinds of studies due to their higher resolution and 
sensitivity compared to ROSAT. We also note that a sim- 
ilar analysis could be carried out in samples of synthetic 
clusters derived from numerical simulations. 
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Cluster 


z 


Exp. time (sec) 


in (kev j 


T ( i n44 „„„ „— 1\ 
Lx(W erg s ) 


A85 


0.0518 


10240 


6.2ot;i:ii; 


19.52+i;^ 


A478 


0.0881 


21969 


6-90lg; 3 3 5 5 


32.00±1S| 


A644 


0.0704 


10246 


6.5912:3 


18.921^:1? 


A1651 


0.0860 


7429 


R in+0.20 
°- lu -0.20 


18.78l|| 


A1689 


0.1810 


13957 


q f)2 +0 ' 40 

3 - ui -() .30 


55.73_ 8 92 


A 1795 


0.0631 


26172 


5.88t° it 


25.42li;« 


A2029 


0.0765 


12550 


8-47l°;« 


41 qQ+ 2 -96 


A2034 


0.1510 


8952 


7.00 


6.86 


A2052 


0.0348 


6211 


o ir>+0.20 
°- ±u -0.20 


4 07+0.34 

-0.34 


A2142 


0.0899 


6186 


9 7D+ 1 ' 30 


61 12+ 3 - 95 
7.09_°;_* 


A2199 


0.0299 


40999 


4 in+o.os 

^■ lu -0.08 


A2219 


0.2250 


11200 


12 40+°' 50 

lz -™-0 .50 


64.561^ 


A2244 


0.0970 


2963 


-0.42 


25.32l»;U 


A2319 


0.0559 


3169 


q 10+0.15 

y - iz -0.15 


39.7411:1? 


A2382 


0.0648 


17444 


2.90 


0.91 


A2390 


0.2310 


10335 


ll.lOlJ;™ 


63.49±1|;|; 


A2589 


0.0416 


7289 


o 70 +1.30 


q /ir>+0.38 
°'^ z -0.38 


A2597 


0.0852 


7163 


4 40+0-40 

^•^ u -0.70 


15.37l^ 9 9 


A2670 


0.0761 


17679 


4 41^+0.20 

^•^°-0.20 


4.97_°*_ 


A2744 


0.3080 


13648 


ii.ool° :» 


62.44li 4 4 :« 


A3266 


0.0594 


13547 


8.ool° :Sg 


16.481°:?* 


A3667 


0.0552 


12560 


7 00+ ' 60 

' - uu -0.60 


22.70l 4 4 : 2 2 ° 


A3921 


0.0960 


11997 


4 qo+°' 55 

^• yu -0.55 


10.92+i'g 


A4059 


0.0460 


5439 


o q7 +0.12 
■J-y_0.12 


5-78l°j 4 


HCG62 


0.0137 


19639 


1 1+0.05 
1 - 1 -Q.05 


0.12 



For A2034 and A2382, the temperatures and X-ray luminosities were taken from Ebeling et al. (1996) 

who do not provide error bars. 
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Table 2. Best fit values for the Sersic parameters of the ICM, k from Eq. (JSJ, and radius r2oo at which the mean 
density is 200 times the critical density of the Universe. 



Cluster v a eq (kpc) n e0 x 10 3 (cm 3 ) k r2oo (kpc) 



A85 


0.55±0.01 


260±5 


5.72±0.11 


7.17±0.32 


2714 


A478 


0.50±0.01 


143±4 


18.51±0.52 


7.04±0.42 


2607 


A644 


0.82±0.01 


389±6 


3.58±0.06 


7.19±0.21 


2421 


A1651 


0.74±0.02 


385±9 


3.57±0.11 


6.59±0.25 


2453 


A1689 


0.58±0.01 


198±9 


14.26±0.82 


7.90±0.35 


2594 


A 1795 


0.54±0.00 


164±2 


12.30±0.23 


7.64±0.21 


2462 


A2029 


0.49±0.01 


151±4 


16.63±0.52 


8.65±0.43 


2946 


A2034 


1.00±0.03 


811±18 


1.72±0.04 


3.58±0.43 


2450 


A2052 


0.47±0.01 


100±5 


11.70±0.72 


9.90±0.67 


1930 


A2142 


0.81±0.01 


528±6 


3.86±0.05 


4.26±0.50 


2824 


A2199 


0.60±0.00 


180±1 


6.82±0.06 


9.42±0.20 


2142 


A2219 


0.85±0.02 


689±15 


3.40±0.08 


4.04±0.24 


2887 


A2244 


0.56±0.02 


237±17 


6.38±0.68 


13.05±0.70 


3039 


A2319 


0.80±0.01 


755±13 


1.86±0.04 


5.13±0.10 


3294 


A2382 


1.17±0.04 


831±19 


0.49±0.01 


6.61±0.66 


1900 


A2390 


0.59±0.02 


313±16 


8.61±0.48 


5.13±0.62 


2658 


A2589 


0.72±0.02 


312±8 


2.31±0.07 


9.25±1.73 


2008 


A2597 


0.39±0.01 


37±4 


71.44±46.76 


12.42±1.60 


2035 


A2670 


0.52±0.03 


215±27 


4.05±0.74 


11.44±0.58 


2331 


A2744 


1.35±0.07 


877±25 


2.26±0.07 


4.82±0.29 


2506 


A3266 


1.18±0.01 


935±7 


1.22±0.01 


5.66±0.20 


2973 


A3667 


0.89±0.01 


990±9 


1.07±0.01 


3.29±0.36 


2723 


A3921 


0.81±0.02 


653±15 


1.41±0.04 


3.66±0.50 


2090 


A4059 


0.64±0.01 


233±8 


4.34±0.18 


8.46±1.02 


2033 


HCG62 


0.36 


22 


17.4 


36.9±1.63 
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Table 3. Gas and dynamical masses, potential energy, central gas column density and integrated specific entropy for 
gas and DM as derived from Sersic parameters. 



Cluster 


M 9as (xlO 14 M ) 


M Dj/n (xlO 15 M ) 


U pot x 10 58 (kg m 2 s~ 2 ) 


So (kg m- 2 ) 


Gas Spec. Entr. 


DM Spec. Entr. 


A85 


5.11±0.41 


2.59±0.35 


0.78±0.20 


0.36±0.01 


204.3±0.18 


206.1±0.91 


A478 


4.98±0.56 


2.29±0.35 


0.72±0.21 


0.68±0.03 


204.1±0.24 


205.7±1.00 


A644 


1.80±0.13 


1.18±0.16 


0.39±0.10 


0.26±0.01 


201.3±0.15 


203.7±0.87 


A1651 


2.53±0.29 


1.46±0.24 


0.46±0.14 


0.27±0.01 


202.3±0.25 


204.4±0.98 


A1689 


4.11±0.85 


2.35±0.52 


1.05±0.41 


0.65±0.05 


203.3±0.43 


205.3±1.22 


A 1795 


3.10±0.20 


1.63±0.20 


0.45±0.10 


0.50±0.01 


203.0±0.14 


204.9±0.88 


A2029 


6.14±0.69 


3.26±0.46 


1.26±0.31 


0.66±0.03 


204.6±0.25 


206.5±1.00 


A2034 


4.44±0.44 


1.79±0.42 


0.64±0.29 


0.23±0.01 


203.1±0.22 


204.6±0.93 


A2052 


1.67±0.38 


0.95±0.21 


0.13±0.05 


0.32±0.02 


202.4±0.49 


204.3±1.31 


A2142 


5.11±0.22 


2.17±0.40 


0.94±0.35 


0.38±0.01 


203.4±0.09 


205.0±0.80 


A2199 


1.27±0.04 


0.86±0.08 


0.18±0.03 


0.28±0.01 


201.2±0.07 


203.5±0.79 


A2219 


8.55±0.85 


3.56±0.75 


2.17±0.89 


0.42±0.01 


204.4±0.22 


205.9±0.94 


A2244 


3.92±1.45 


3.33±1.15 


1.53±0.88 


0.36±0.05 


203.7±0.75 


206.4±1.60 


A2319 


7.45±0.62 


3.64±0.63 


1.79±0.61 


0.26±0.01 


204.5±0.19 


206.4±0.90 


A2382 


0.96±0.01 


0.67±0.11 


0.11±0.03 


0.06±0.01 


200.4±0.22 


202.9±0.91 


A2390 


9.06±2.17 


3.72±1.01 


1.78±0.85 


0.62±0.05 


205.0±0.51 


206.5±1.31 


A2589 


0.97±0.12 


0.73±0.11 


0.13±0.04 


0.15±0.01 


200.6±0.28 


203.2±1.00 


A2597 


2.30±2.14 


1.26±1.13 


0.22±0.49 


0.78±0.52 


203.1±1.69 


204.9±1.94 


A2670 


2.86±2.09 


2.03±1.35 


0.46±0.59 


0.22±0.05 


203.6±1.24 


206.0±2.18 


A2744 


4.00±0.52 


2.20±0.47 


1.36±0.55 


0.26±0.01 


202.3±0.28 


204.4±0.95 


A3266 


3.33±0.11 


2.03±0.31 


0.92±0.28 


0.17±0.01 


202.4±0.07 


204.7±0.77 


A3667 


6.93±0.31 


2.52±0.57 


0.83±0.38 


0.19±0.01 


204.4±0.10 


205.7±0.81 


A3921 


3.53±0.40 


1.34±0.31 


0.29±0.13 


0.17±0.01 


203.1±0.25 


204.5±0.98 


A4059 


1.24±0.20 


0.80±0.14 


0.15±0.05 


0.22±0.01 


201.2±0.34 


203.5±1.09 



